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^Q ' Abstract. An invariant Gibbs' state for the nonlinear Schrodinger equation on 

pH I the circle was constructed by Bourgain, [B] , and McKean, [MC] , out of the ba- 

sic Hamiltonian using a trigonometric cut-off. The cubic nonlinear Schrodinger 
equation is a completely integrable system having an infinite number of addi- 
tional integrals of motion. In this paper we construct the second invariant Gibbs' 
state from one of these additional integrals for the cubic NLS on the circle. This 
^ ' additional Gibbs' state is singular with respect to the Gibbs' state previously con- 

0^ ' structed from the basic Hamiltonian. Our approach employs the Ablowitz-Ladik 

^ . system, a completely integrable discretization of the cubic Schr5dinger equation. 

a^ 
o 
o 
o. 

1. Introduction. The nonlinear Schrodinger equation for the complex function 

-0(0;, t), X G S, t G -R"'^ is 



i^« ^_^"+p|^|2p-2^^ 



> 

X 

j^ ■ where p is an arbitrary integer. It can be written in Hamiltonian form. Let 

A^ = ('05 "0) be a space of pairs of two complex functions ifj = Q + iP and 
t/; = Q — iP on the circle § of perimeter 1. For any two functionals F and G on 
A4, define 
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2 K.L. VANINSKY 

The Hamiltonian 7i{'i(j , t(j) = J^ |V''P + li^l'^^dx produces the nonhnear Schro - 
dinger flow 

The equation has two other generic integrals. They are A/" = / \'i(j\'^dx, the 
number of particles, and V = — J it/j't/jcix, the momentum. 
The Gibbs' measure associated with the basic Hamiltonian 

1 1 

e~^'^dvol = e~2'^ — - Alv = e~^'^ — r /\ i6ip(x) A 6'iIj(x) 
oo! ' ' oo! ' ' 

oo x6S 

was constructed by Bourgain and McKean [Bl-2, MC]; it is the product of two 
independent copies of circular Brownian motion for the components of the func- 
tion ijj = Q + iP coupled together by the nonlinear factor e~ 2 Js '^ +^ I '^^. This 
measure is invariant under the flow; for any p the later exists almost everywhere 
with respect to the Gibbs' measure. 

The cubic NLS corresponds to the case p = 2. It has inflnite series of com- 
muting Hamiltonians. The flrst flve are listed below 



Hi = / \ij\^dx, 
H2 = / —iip'ipdx. 



The classical generic integrals A/", V and Ti are the flrst three of this inflnite 
series. All of them are integrals of the isobaric polynomials, counting if) and d/dx 
as having degree 1; for example, Ti^ is an integral of an isobaric polynomial of 
degree 6. 

We consider the phase space M. of pairs (^, ^) from the Sobolev space H^* . 
The original equation deflnes a global flow on A1, [B2] . In this paper we construct 
an invariant Gibbs' state from the higher Hamiltonian Ti^: 



e ^'^^dvol = e a'^s — - A i5ijj(x) A 5ip{x 



00! 



If ^{x) = J2k e^ '>P{k), then the Sobolev norm of ^{x) is \xP\jj, = J^k l^('^)r(l + ^ )' 
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This state is singular with respect to the Gibbs' state e~'2'^^dvol. To construct 

the measure and prove its invariance we wiU use the Ablowitz-Ladik (AL) system. 

The AL equation for the complex N-periodic function i/j{n,t) ,n E Z^, t E R^, 

is* 

It is easy to check that the quantity 

N-l N-1 



D=i[Rk= n(i-iV'fci 



A:=0 k=0 

is an integral of motion. We assume that \t(jn\ < 1 for all n. This area of the 
phase space we call the box B; it is invariant under the flow. In the box all 
quantities Rn and D are positive. 

In fact, AL has many integrals of motion. The first three interesting ones are 
-^0, h, Ia:'- 

N-l 



Nl2= Y.'^r^r-l. 



r=0 

^-^ _ I _ 

r=0 

From these integrals we form Hamiltonians 

ifi = N{Io + /o) 

Hs = N{h + h- 2/o - 2/o) 

ifs = N{h + I4- 4/2 - 4/2 + 6/0 + 6/0). 

Let Mtv be a space of complex A^-periodic sequences ifjn = Qn + iPn, n E Z^ 
with li/jjil < 1. Introducing the bracket 

^ ' Sf Sg 6f Sg 






6ijn 5^n 5'^r, 



for two functionals / and g on M;^ we can write the original AL flow in Hamil- 
tonian form 

•^n = {'i/'n,i^3}^„- 



denotes a time derivative; ijjn = V'('^i Oi periodicity means xpn+N = '4'n for n £ Z. 
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The symplectic form loq producies the volume element 



The volume of the box is infinite 



dvol= ,^, ^ 

NIDn- 



dvol = I „ /\ iSipn A 6iJn = oo 



NIDn 

B B ^ 

due to the singularity of the factor jy—. On the box we can define the Gibbs' 
state with the density 

_ jv£ TT _ jv£. rr 1 A _m1. U 1 Ac-/ f T 

e 2^5^vol = e 2^^^—p^u)Q=e ' n\D /V^'^^^^'^n, 

n ' n 

The singularity of the volume element near the boundary is rectified by the 
vanishing factor e-^-isJVJo _ ]j3N . ^j^^ total mass 

_ JV£ TT 

e 2 -"sdvol < oo. 

B 

The area of the box close to the boundary has negligible probability. The measure 
of the area outside the box we take to be zero. 

Let e = jj. Under the scaling ipn = ^i^{j^) and 9 = e~^t, the AL vector field 
approaches the NLS vector field as A^ -^ oo. Indeed, 

V'n+l = ^n + ^e + 2^n^^ • • • . V'n-l = ^l^n ' ^l^'n^ + -ll^n^"^ ■■■■ 

So, 







V'n+l + i^n- 


-1 - 


-2V'n = 


-i^'y 


+ .... 


and the fiow in the time scale 6, 












' 89 = 


-- -(V'n+l +V'n- 


-1 


-2V'n) 


+ |V'n 


l'(V'n+ 


scales 


as follows: 













ot 

This indicates that AL flow converges to the NLS flow, though statement can 
not be taken literaly because these two flows live in two different spaces Mat and 
M. 
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The integrals Hi, H^ and H^ of the AL system converge to Tii, H3 and H^: 

iVifi(V'n,V'n)-^Hl(^,^), 

Also, for any i/j, we have DN^i/jmi^n) -^ 1, as A^ ^ 00. Therefore, one should 
expect to get from the Gibbs' state for AL, in the limit of A^ ^ 00, the desired 
invariant measure 



^'^''dvol^e a'^^ — - /\ i5ijj{x) A 5ijj{x) 



e 2' 

00! 

Again, this statement is correct but can not be taken literaly, because these 
measures do not live on the same space. To make the arguments rigorous we 
need to embed the AL flows and measures into the function space. For this, we 
use interpolating trigonometrical polynomials. 

The main goal of the paper is the space-time random field -(/'(a;, t), a; G §, t G 
B} such that: 

(i) -(/'(x, t) is stationary respect x and t; 

(ii) ip{»,t) has Gibbs' distribution e~2'^^dvol] 

(iii) the random variable il){»,t), t 7^ is measurable with respect to the 
cr-field generated by '^{•^ 0); the measure is supported by the solutions of NLS; 

(iv) The x-derivative of almost every realisation of the random field is Holder 
continuous: 



\dj,ip{xi,ti) -9xV'(a^2,t2)| < K 



Fl - iC2 P + \tl - t2h 



with a random constant K. The exponents ^ and ^ are optimal. 

The paper is organised as follows. In the section 2 we introduce the commut- 
ing fiows of the AL hierarchy and commuatator formulas for them. Section 3 is 
devoted to the study of the direct spectral problem for the auxiliary linear sys- 
tem. We explain the details of the spectral curve for various types of potentials. 
Invariant quantities are computed in section 4. The Floquet and dual Floquet 
solutions are studied in sections 5 and 6. Section 7 provides formulas for the 
symplectic structure and the Poisson bracket. This completes the first part of 
the paper about integrability properties and Hamiltonian formalism for the AL 
system. 

The second part starts with section 8 where the strategy for constructing the 
measure is outlined; it is implemented in the subsequent sections 9-11. 

Finally I would like to thank I. Krichever, H. McKean, V. Peller and J. Zubelli 
for helpful discussion. It is also pleasure to thank MPI in Bonn and IMPA for 
their hospitality. 
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2. Ablowitz— Ladik hierarchy. The Ablowitz-Ladik equation for the complex 

A^-periodic function "ipin.t) n E Z^ , t E R^ is 



(1) ii^r 

Recall that 



-(•071 + 1 +•0^ 



n-1 



2^n) + |'i/'n|^('i/'n+l + V'' 



n-1 



7V-1 JV-l 

A;=0 fc=0 

is an integral of motion, whence the "box" B — [t/; : {t/jnl < 1 for all n) is 
invariant under the flow. In the box R^ and D are positive. 

The flow (1) is one of infinitely many flows of the AL hierarchy. The rotation 
or the phase flow 

• 

(2) iV'n = — 0n 

is the flrst flow of the hierarchy. Let 



V2{n,t,X) 



1 



'Rr 






The phase flow is the compatibility condition for 

(3) [dt-Vi,A-V2]=0, 

where Vi = i(Ts/2 and A/^ = fn+i is a shift operator. Here and below a denotes 
the Pauli matricies 



0-1 



1 

1 



0-2 



-i 
i 



0-3 



1 
-1^ 

The original AL equation (1) is the compatibility condition for* 

(4) [dt-V3,A-V2]=0, 

where 
Vs{n,t,X) 

A^ - 1 + I'i/'n'^n-l - iV'nV'n-l 'i/'nA - -^n-l A"^ 

•i/'n-l-^-V'n-^"^ 1 - A-2 + i'i/'nV'n-l + l^r. 

The formula means that** 

{dt - Vs{n + 1))(A - V2{n)) - (A - V2{n)){dt - Vs{n)) = 

i.e. 

V2{n) = Vs{n + l)V2in) - V2{n)Vs{n). 



*A similar form of the commutator formalism was considered in [AL, MEKL]. Our form 
has the small advantage, that it leads to a unimodular monodromy matrix. 



^The operator V = V2 acts like shift: V2(n)f, 



n — '-n 



-I- 
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3. Monodromy matrix. The spectral curve. Since i(^n+N = i^m all ma- 

tricies Vfc, /c = 1, 2, . . ., satisfy the periodicity condition Vk{n + N) = Vk{n). We 
introduce the transition matrix 

Tn,^{t, A) = V{n - 1, t, X)V{n - 2, t, A) . . . V{m, t, A) {n>m). 

It is easy to see that the spectrum of Tm+N,m{tj X) does not depend on m. 
Indeed, 



-m+l+N,m+l 



(t, X) = V{m + N)... V{m + 1) = V{m + N)Tm+N,mV{m)-\ 



Also, the spectrum of T/v = Tn,q does not depend on time. This follows from 
the identity T^ = ["1/3(0), T/v], which is proved as follows: 



dt[V...V{k + l)V{k)...V] 

= ... + V...V{k + l)V{k) ...V + V...V{k + l)V{k) ...V + ... 
= ... + V...{V3{k + 2)V{k + 1) - V{k + l)Vs{k + l))V{k) ...V 

+ V...V{k + l){Vs{k + l)V{k) - V{k)Vs{k)) ...V + ... 
= Vs{N)Tn{X)-Tn{X)Vs{0). 

Consider the special "Floquet" solution 



9nW 

9lW 



of the eigenvalue problem g„_|_i(A) = V{n, A)gn(A) specified by the condition 

gAr(A) = TArgo(A) = tygo(A), 

where w being the complex "multiplier" determined by 



= det 



rp21 



rpl2 

TW-w 



zi;2 _ 2'u;A(A) + 1, 



where A(A) = trT/v(A)/2. We know that A(A) is an integral of motion. The 
multiplier is 

w = A(A) + VA2(A) - 1; 

It becomes single-valued on the spectral curve 
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The curve F inherits symmetries from V2{n, A). It is easy to check that 

1 



(1) a^V^=ja^ = V{X), 

(2) crsV{-X)a3 = -V{X). 
Obviously T/v(A) satisfies similar identities. Whence 

1 



(3) A ^=j = A(A), 

(4) A(-A) = (-l)^A(A). 

The symmetries (3-4) produce an antiholomorhic involution r^ and a holomor- 
phic involution r on the curve 

(5) Ta : {X,y) — ' [yy 

(6) r: (A,y)^(-A,(-l)^y). 

On r, there exists another holomorphic involution t± which permutes the sheets, 
viz. T± : (A,y) — s> (A, — y) . From the quadratic equation for w, we have 
w{Q)wiT±Q) = 1. 

Remark. Obviously, we have freedom in the choice of sign for the second 
coordinate of the involution, say 

Ta : (A, 2/) ^ ( T'^^ 

The sign in (5-6) is chosen in such way that r^ and r preserve the infinities of 
the curve (see example 5 of this section). 
Example 1: vanishing potential t/jn = 0. Then 

Tn{X) = (^^ ^^N^ , A(A) = coshATlogA and w{Q) = e±^i°g^(«). 

The points A^ = e'a^'^, /c = 0, • ■ ■ , 2N — 1 satisfy the condition A^(A^) = 1 and 
at these points* A*(A^ ) = 0. The points A^ are simple crossings of the curve 
F. They form periodic/antiperiodic spectrum; namely, w = +1 for k even and 
w = -1 for k odd. If A = e*^ then A{e'^) = {e'^^ + e"*^^) /2 = cos^iV. The 
graph of A(A) for A^ = 2 is shown in fig 1. 

When 'i/'n is not identicaly zero, then the double points A^ split into pairs A^^ 
and A^ of simple ramification. The symmetry (3) implies that A(A) is always 
real for |A| = 1. Note also that, due to the symmetry (4), the branch points form 
symmetric pairs A^, A^ and i Aj^,, A^, such that Xj^ = —Xj^,, when k — k' = 
(mod N). 

The branch points Aj^ lie only on the unit circle**. Indeed, the equation 



*• now denotes derivative in A variable. 

**The statement is proved by the method of [AL]. 
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FIG. 1 
^_l_l(A) = F(n, A)gn(A) can be written in the form A^g„(A) = Ag„(A), where 



A. 



-A 



VRnA 



and Afn = fn+i- 



It is easy to see that A is unitary in the space of 2A^-periodic vector functions with 
the complex inner product < f , g >2<p= j}^J2k=o fk9k- Indeed, introducing 
the formal inverse 



A: 






it is easy to derive the Cauchy identity 

< Af,g>c? =<f,A-^g>(p + 

(7) 



+ -Tf^/^N-l {fN9N-l - fo9-l) + -TtV^N-I {f-i9o - fN-l9N) 



N 



N 



where < f , g >(?= -^ X]fc=o fk9k- Appling this formula twice, first to the interval 
/c = 0, ■ ■ ■ , A^ — 1 and then to the interval k = N, ■ ■ ■ , 2N — 1 we obtain a similar 
identity for < •, • >2^ which implies the result. 

Example 2: A^ = 2-periodic case. It is instructive to analyze this case com- 
plitely. The potential can be written in the form 



V'r 



A + B, n is even 
A — B, n is odd. 



10 
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The matrix T2 can be easily computed: 



X^+{A-B){A^B^ X{A + B) + X-\A -B) 
X(A-B)+X-\A + B) X-^ + {A + B)(A-B) 



ThenA(^) = ^ [cos2^ + |A|2 - \B\'^] for A = e*^ and D = {1-\A+B\'^) x(l- 
\A-B\^). Consider the case S = and A ^ 0. Then A(^) = ^ [cos 2^ + \A\^] , 
and D = (1 — |Ap)^. The graph of A(6') is shown in fig. 2. 




FIG. 2 



The double points Aq and A2 split into pairs, while A^^ and A3 remain double. 
Now consider the case when A = and S 7^ 0. Then, A(^) = -^[cos2^ 

— |i?p], D = (1 — l-Bp)^. In this case, the double points Aq and A2 do not split, 
while X^ and A3 split into pairs of simple roots. 

Finally, there is an open area in the space of parameters, say, A ~ S, where 



A(^) 



cos 2^, 



D^ 1-4\A\^ < 1. 



In this case all double roots split into pairs of simple roots. 
Example 3: constant potential V'n = V'o- Consider, first, the case A^ 



2. Then 



To 



1 



A^ + I^oll 
Xipo + A" Vo 



Xifjo + X Vo 

A-2 + l^oP 
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A{e) = -^ [cos 26'+ l-i/'oP] , and L* = (1 - |^oP)^- The branch points AJ are 
given by the equation A^(A) = 1. Obviously, the points A^^ and A3 do not spht. 
For Aq we have the equation 

l-|Vo|' = cos2^+|^o|', 

so for small |V'o|, 

l-2\ijo\^ = 1-29^ + 0(9^). 

Finally, we arrive at 9 = IV'ol + Od'i/'ol^) and 

;^± ^ ^±i{\i,o\+0{\^o\^)) y^± ^ ^±i{\i,o\+0{\^o\^))+i7r _ 

In words, the open gap is proportional to the absolute value of the potential |?/'o|. 
The case of general A^ can be treated easily. One has to compute spectrum 
of T/v. After simple algebra 

x± _ ±*(|Vo|+0(|^o|')) x± _ ^±t{\^o\+0{\^o\^))+i7r 

Aq — c , AjY — e 

Example 4: two-gap trigonometric potential. Consider the A^-periodic poten- 
tial ip'^: 

V'' = e^^'=>o, <^^ = 1^ fc = 0,---,iV-l. 

We will show that this potential opens /c-th and k + A^-th gap as in 

(8) _X^ = e*'^'=/2±^(IV'o|+o(|^o|3))^ 

and 

(9) ;^± = gi<^fe/2±J(|Vo|+0(|V'o|^))+i^_ 

To emphasise the dependence on -0'^ we write V{n, X\i{j^) = V'^{n), T(A|'i/'^) = 
T^(A) and A(A|V''') = A'=(A). Note, first, that 

L''{n) = $"y'^(0)$-'', where $ = e^'"^'^'^/^. 

Therefore, 

T(A|e^'^'=>o) =V''{N-1)--- V''{0) = $^($- V=(0))^ 



12 
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Taking the trace we obtain 

This and the result of Example 3 produce the formula (8). (9) follows from this 

and the symmetry (4). 

Remark 1. The situation here is similar to the NLS equation, [MCVl], where 

the A^-th gap opens in proportion to the A^-th Fourier coefficient. 

Remark 2. The AN branch points Aq , ■ ■ ■ , X2N-1 determine the curve. The 

symmetry (4) leaves only 2N free parameters. In fact, there are only A^ free 

independent parameters due to periodicity conditions. 

Example 5: generic periodic potential. The potential can be written as a sum 

of harmonics i/j'^, /c = 0,---,A^ — 1. All gaps are open. To construct the curve 

r, let us take two copies "+" and "-" of CP^ cut along circular arcs connecting 

X~ and A+ (fig. 3) 





FIG. 3 

Each copy CP^ has two marked points PJ} and Pr, . The behavior of 



0/00 



O/oo- 



A(A), y{Q) = a/A(A)^ — 1, and w{Q) near these points is this: 

A(A) ~ ^A^, y{Q) ~ ±^A^, w{Q) ~ A±^, A = A(g), Q e (P±); 



and 
A(A) ~ ^A-^, yiQ) ~ ±^A-^, wiQ) ~ A^^, 



A = A(g), Qg(Po^). 



From this, we see that r^ maps point Q = (A, y) in the vicinity of Poo to the 
point TaQ = {k,y) in the vicinity of Pq . The curve is obtained by gluing 
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13 



P^+ P^= P^ 




real ovals 



P +P =P 

00 



FIG. 4 



together "±" copies of CP^ along the cuts and identifying the points Poo with 
Pq as specified 

P+=P++P+, and P-=P-+P0"- 

The cuts are the projections on the A-plane of the real ovals of anti-holo- 
morphic involution Ta- The infinities P"*" and P~ are also fixed points of Ta- The 
multiplier w{Q) is single- valued on F and satisfies the identity w{TaQ) = w{Q). 
This guaranties continuity of w{Q) at the points P+/~, where w{Q) — > oo/O 
respectively. The involution r is defined so that it preserves the infinities. The 
formula w{Q) = A{X{Q)) + y{Q) implies w{tQ) = {-l)^w{Q). 

Let us introduce a multivalued function p{Q) on F by the formula w{Q) = 
qp{Q)n _ Obviously, p{Q) is defined up to the integer multiple of 2'iii/N . The 



14 
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differential dp is of the third kind with poles at P^ and P~ . The relation 
w{Q)w{r±Q) = 1 implies that p{Q) = —p{t±Q) (mod 2Tvi/N). Therefore, in the 
vicinity of P^ 






p+(A) = +logA + /o + y + y + ■■■ , 
P^W = -logA-/o- y - y H , 

where A = X{Q). The fact that r preserves Poo and w{tQ) = {—l)^w{Q) 
implies Ik — for k odd, so 



p+ (A) = +log A + /o + -^ + -T + 



Poo(-^) = -logA-/o 



A^ 



A2 



A4 



A4 



+ 



The anti-involution Ta maps Poo to Pq . This and relation w{TaQ) = w{Q) 
imply, in the vicinity of Pq , 

p+(A) = -logA + 7o+72A2+74A4 + ---, 
p^(A) = +logA-7o-/2A2-74A4 + ---. 

The explicit form of the integrals /q, /2, ■ ■ • will be computed in the next section. 

4. Integrals of motion. To simplify the calculations, we introduce Vnew{n) — 
^fPnYoiAip)- The corresponding mulitipliers are Wnew = V^tCoid, and 

Pold = Pnew - — log P*. 

We write all formulas in this section for the new spectral problem. 

We introduce T^,^ = (l+Wn(A))e^"'"^(^)(l+Vr,„(A))"^ where Z is a diagonal 
matrix and W is antidiagonal. Separating the diagonal and antidiagonal part in 
the equation Tn+i,rn — V{n)Tn^rm we obtain 



(1) 

(2) 



W. 



„+i,^ ^"."^ = A + \1/„1^, 



+ie 



AW^n + ^, 



where 



A 



A 
A-i 



and ^, 



_0 Vn 
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Let Wn{X) = E ^nA"^ where 
k=i 



w: 



w'^+ 
w'^~ 



From (1-2), we have Wn+i{A + ^nW^n) = AW^n + ^n- Introducing 

and (T_ 



(j+ = 



1 





1 



we can write A = cr_|_A + (T_A ^ and 

J]vr^\iA-'=(Aa+ + A-V_)+ J] W;+iVi>^W^^PA- 



■r—p 



(3) 



fe>i 



A;,p>l 



fc>i 



Collecting terms with A^ in (3), we have VK^_,_]^a_|_ = cr+VF^ + \E'^ and 

Collecting terms with A~^ in (3), we have W^_^_ia+ = a+W^. This implies 
W^ = 0. 

Collecting terms in (3) of the order X ^, k > 2, we have the reccurence 
formula 

(4) w!:+la+-a+W::+'=a-Wt'-W::-la-- J^ VT^+iVEt^VT^P; 



r+p=k, 
r,p>l 



whence, for /c = 2, 



W^^ 



-t/jn+l +ll^n+l\A 



V'n-2 -'0n-2lV'n-lP 



Therefore, for n = N and to = we have 



Tn{X) = (1 + W^(A))e^--o(^)(l + Wo{X))-^ 



and 



AT-l 



e^~."(A) = JJ(A + ^^VF'^ 



r=0 



16 
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Using the fact that \E'rW^r is diagonal we infer that for Q G [P^] 

N-l / 



r=0 



k>l 



Therefore, 

p+(A)iV = iVlogA+ 



^ N-l _^ /N-l .^ 



r=0 \r=0 

Using the exphcit expression of W^, W^, we have* 

N-l 

NI2 = Y, i^ri^r-l, 
r=0 

^-^ _ 1 _ 

r=0 

5. Floquet solution. Baker— Akhiezer function. The Floquet solution was 
defined in section 3 as a special solution of the spectral problem gn+i = V{n)gn 
with the property gN = ^Vgo = wgQ. It is uniqlue specified by the boundary 
conditions g^^ + gl\n=o = I- 
Example: vanishing potential t/jn = 0. Then, 

A 



V^2(n,A) 



A- 



and Tn{X) 



A^ 
A-^ 



Therefore,** for Q e {P+/P^), 

gn(Q) = e±-^°^^[go+/goT 
andforQG(p-/Po-) 

gn(Q) = eT-^°^^[go-/go-] 



«;(Q) = e±^i°g^; 



«;(g)=eT^^°sA. 



*We write the coefficients for the old spectral problem. 



•"For any g 



we write g 
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go 5 


go are 














+ 


1 











go = 





go = 


1 



We introduce another Floquet solution g^ normalised by the condition i 






g;. 



i^Z |n=o = 1- Let us define 



n-l 



Do = l, Dr,= \[Rk, n = l,...,iV-l; Dn=D. 



k=0 



Lemma 1. i. The Floquet solution gn{Q) is: 



with 



.{Q) = C{Q) 



C{Q) 






+ {l-C{Q)) 






A = A(g). 



^2 



w{Q)-T 



'22 



T12 _ Til + «;(g) T2i-T22 + «;(Q)^ 

here, T = T/v . 

a. The Floquet solution gn{Q) is 



with 



kn{Q)=R{Q) 



R{Q) 









iT 



12 



+ {i + R{Q)) 

iT^^ - iw{Q) 



w{Q)-T^^-T^^ w{Q)-T^^-T^^' 
Hi. The following formulas hold: 



and 



CTlgn{TaQ) = gn{Q), <yign{Q){raQ) = gn{Q)- 



gn(TQ) = {-iriasgniQ). 



(iv). The Floquet solution gn(<3) has 2N poles at the points 7's on the real 
ovals of the curve T. On each oval there is just one 7. 
(v). If Q E {P^/Pq), then gn(<5) has the development 



with 







-1 


00 


00 


1 


inlogA -L 


5^g+HA-7$^g+(n)A^ 


y ^n 


Ls=0 s=0 J 


+ 
— 


"1" 




1 


g^(n) = 


■ -V'-I " 





18 



K.L. VANINSKY 



Also for Qe{P-/Po), 

gn(Q) = e^"^°s^/D; 
with 



Y,S7in)X-yJ2^7in)y 



.8 = 



h (n) 



s=0 

^0 



Proof, i. Obviously, 



ln{Q) = C'iQ) 



" n 



TlH\)i 



+ C"{Q) 






with some C'{Q) and C"{Q). The normahsation (7^ + gl 



\n=0 



1 imphes 



C"{Q) = 1 — C'{Q). Another condition T„go = wgo produces the system 



rpll rpl2 

-'-N -'-N 

rp21 rp22 

-'-N -'-N 



c 

1-C 



w 



c 

1-C 



Solving the system, we obtain the stated formulas for C{Q). 

ii. The proof is identical to the proof of (i). 

iii. The proof is based on explicit formulas of (i) and (ii). Using (1) of section 
(3), we have 



Tn,o(A(raQ)) = (TiT,,o(A(Q))ai 



and 



w{TaQ) = W{Q). 



From this we have 1 — C{TaQ) = C{Q). The proof is finished by substituting 
Tni^TaQ)) and C{TaQ) into explicit formula of (i). Similarly R{TaQ) = i+R{Q). 
This implies the stated formula for g„. 

The last formula for gni'^'Q) is proved along the same lines using (2) and (6) 
of section 3. 

(iv). It is easy to see, by perturbation arguments, that fXk = ^(7/0) are near 
the points A^ = e^^^ ^ /c = 0, . . . , 2N — 1 for a small potential. We will prove 
that |/Ufc| = 1 always. 

Let gn(Q) have a pole at some point 7. Then, from the formula of (i), we 
have 

T^^-T^^ = w{Q) and T^^ - T^^ = w{Q). 

Consider the special solution fn(A*) = T^ (//) — T^ '(n) of the eigenvalue problem 
Anfn = A*fn- We have 



1 



and 



In 



w{Q) ■ 
-w{Q) 
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Now the Cauchy formula (7) of section 3 to obtain 

< Af , f >(p=< f , A~^f >cp +boundary terms. 

Note, that Af = ^uf, A~^f = n~^f. It is easy to compute the boundary terms: 
to wit, 

l(|«;|2-l)(^-i-^). 
Substituting this into the Cauchy formula, we obtain 

1 



(1-i^r) 



A^ 



kl -l)-(f,f) 



0. 



This implies \n\ = 1. 

(v). Consider P^/Pq*". The relation between two asymptotic expansions 
follows from the formula aign{TaQ) = gn(Q)- The actual form of the coefficient 
is derived by substituting the asymptotic expansion 



.(Q) 






^Dr 



s=0 



into the spectral problem g^j,-!-! = V{n, Xjgn, where F(n, A) = [Aa_|_ + A ^cr_ 
+\1/„1. We arrive at the reccurence relation 



The boundary condition (7^ + Qnln^o 



1 implies 

^o('^) + ^o(^)U=o = 1; glin) + 9^{n)\n=o = for s > 1. 

Starting from s = —1, reccurently, one can compute gQ"(?i), gf{n), etc. D 

On the curve F one can consider the Baker- Akhiezer (BA) function g(r, n, t. A) 
with 2N poles on the real ovals. On each oval there is just one pole. The function 
g(r, n, t. A) has the following assymptotics at the infinities: 



62 



g(r, ?i,t. A) 
g(r,?i,t. A) = e~i 
g(r,n,t. A) = e~i 



T+n logA+i(A^-l)t_ 



-. 00 



s=0 



aogA-i(A"^-l)t 



-nlogA-i(A^^-l)t 



/d^Es^"^" 



s=0 

CXD 






CXD 



g 



(r,n,t,A) = e^+->°s^+^(^^-^)*V^5^grA^ 



Qe(P+), 

Q e (Po+), 
Qe(^o")- 



s=0 
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The BA function with these properties exists and defined uniqly. The BA func- 
tion g(T, n, t, Q) satisfies the identities 

[dr-Vi]g{r,n,t,Q) = 0, 
[A-y2]g(r,n,t,Q) = 0, 
[at-V3]g(r,n,t,g) = 0. 

The exphcit form of the V^s is given in section 2. The Floquet solution is a 
particular case of BA function with the variables r, t fixed. 

6. Dual Floquet solution. Variational identity. One can write the equa- 
tion V{n)gn = gn+i in the form* 



[J A - JV{n)]gr, = 0. 
Let us define the dual Floquet solution g+ = [(7^^, (7^^] at the point Q by 

g+(Q)=gn(r±Qr. 
Lemma 2. The dual Floquet solution g^iQ) satisfies** 

-{Q)[JA-JV{n-l,X)] =0. 



r+l 
3n 



Proof. The equation l/(n — l)g^_i 
where 



V-\n-l) 



[n can be written as gn-i = V ^(n — l)j 



or, in coordinates, 



9n-l 



9I-I 



1 



Rn-1 
1 



V -Rn-1 

We can also rewrite it in the form 



i9n^ 9n) 



A 
-A 



(A ^gi-tl)n-i9l), 
{-^n-i9l + ^9l)- 



i^n-l A ^ 

-A -V'n-l 



Rn-1 



0. 



This is exactly the stated identity. D 

We need the standard formula [KP] for variations of the quasi-momentum. 



*J = ia2. 

**f A — f , 
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Lemma 3. i. The expression g^(<5) Jgn(Q) does not depend on n and 

g+(g)jgn(Q)=<g+jg>=vi/(g). 

a. The following identity holds 

6p^{Q) =< g+J6V{n - l)g^_i > . 

Proof, i. Can be checked using difference equation for gn(<5) and g^(Q). 

ii. Denote by V{n,X) and gn(A) deformed matrix V{n,X) and the Floquet 
solution gn(A); to wit, V = V + e6V + o(e) and g = g + eSg + o(e). Then, 

g+(g) [(JA - JV{n - 1, A))g„_i(Q)] = 0, 

[g+(Q)( JA - JV{n - 1, A))] gn-i(Q) = 0. 
Subtracting, we obtain 

iV-l 

^g+(JAg._i)-(g+JA)g._i 

n=0 

N-1 

= J] g+ (jy(n - l)g^_i) - {g+JV{n - 1)) g^_i. 

n=0 

It is easy to see, that 

N-l 

RHS = e 5^ g+JdV{n - l)gn-i + o{e). 

n=0 

Using the formula 

7V-1 N-l 

Y, (g+ JAg)g,_i = Y, Si(J^Sn-l) + g+i Jg-1 - g^_i JgN-1, 
n=0 n=0 

we obtain 

LHS = g^_i JgAT-i - g+i Jg-i. 

From the definition of the Floquet solution 

g^_i=e-^Pg+i and gA._i = e^^'g-i, 
we have 

LHS = (^e^(p-p) - l) gli Jg_i = eN^pgt^Jg-l + o(e). 
Collecting terems with e, we obtain the stated identity D 

Let us introduce g* (Q), the dual Floquet solution normalized by the condition 
< g* Jg >= 1. Obviously, 

g* (Q) has poles at the branch points of the curve. 
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7. Hamiltonian formalism for the Ablowitz-Ladik system. As shown in 
[KP], the formula 



defines a closed, nondegenerate 2 form on the space of operators A — V2 with 
periodic potential. Our goal to compute the residue explicitly. 
Near P^, 

(1) ^ = ^° + ^i^' + '""' Qe(Po+). 

The identity 

^(r±Q) =< g+(r±Q)Jg(r±g) > =< g^(Q)Jg+(Q)^ > 

= - < g+(Q) Jg(Q) >= -^(Q) 
implies 

(2) ^ = -^o-^i^'----' Q^iPo)- 

The involution r^ maps Poo to Pq , and gn{TaQ) = o"ign(Q)- Therefore, 
using r±ra = TaT±, we have 

^(TaQ) = g^{r±TaQ)Jg{raQ) = g^ {TaT±Q)Jg{TaQ) 

= g^{r±Q)alJa,g{Q) = -^0) 
so that, if Q e (P+), then TaQ G (Pq^) and 



(3) ^(Q) "^{TaQ) y^^+Xl^'") 

= -V'o-|^----, Q^{P^)- 

Again, applying r±, we have 

(4) ^ = ^,„ + |i + .... Qe(P-). 



AN ADDITIONAL GIBBS' STATE 23 

It is easy to compute 



ijo =< go JSo >= 1 

ipi =< ko^Jgt > + < 

Similarly, introducing S{Q) =< g^{Q)JdV{n — 1) A dgn-i{Q) >, we obtain 



^1 =< go ^Jgt > + < gi ^"^go" >= ^0 - ^-1, etc. 



(5) S{TaQ) = -S{Q). 

One can show that, in the vicinity of the infinities: 



S{Q)=co + ciX + ---, Qe{P+) 

'X 



S{Q)=do + ^ + ---, Qe{P+); 



and 



S{Q) = ao + aiX + --- Qe{Po), 

S{Q) = bo + ^-f + --- Qe{p-). 

(5) implies d-k = —ct and b-k = —dk- Computing residues, we have 

b dX . \ / T \ dX J 

res — -- = res (do H ) (-^o ) -^ = "oV'o, 

p+ y/ A p+ A 

b dX . \ / 7 \ dX J 

res — — = res (6o H ) [Wo ^ j ^ = -OoV'o, 



^ A p£ ^ "^ ' ^"^ ^ A 



and 



Similarly, 



Eb dX 



and 



P+/- 

■'■ no 



res — — = res (cq H ) (V'o H ) -;- = CqV'o, 

P+ ^f A p+ A 

b ctX , \ / I \ dX 

res — -- = res (ao H ) {-^jo ) -^ = -aoV'o, 



Xl ^®^ ^X ^ ^° ~ "° ^ ~^^^ ~ ^^'^ 



^0 



+ /- 
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Using the formulas for Floquet solutions from Lemma 1, we have 

+T 



+ < 



'n 

'0 -.J 5 '"-' A5JDr 



Dn, A/-Rn-1 



\1/ 1 

^ n — 1 



Igl > 



.+T 



^+ 



+ < ^J6^^ A 5v/A^gr > 



+ < 



Dn ^JRn-l 

'1 ,J5 "ll^Z^ AJV^ 



-Dn V-Rn-1 



-iSo > 



2 2 

< ^^5^^^ A 6JD^_i > 



Dn jRn-\ 



+ < ^(5 J A 5v/A^^n-l > 



- < 



V-Dn \jRn-l 



After simple algebra, 

1 ^ SRn-l . 5Dn-l 



.'.= T < 



A 



A5Vl^n-l > • 



1 



Similarly, 



. ~ D ■ • n > - < ^ (^^n-l A 5ijri-l > ■ 

4 -Kn-1 -L^n-l -Kn-1 



rfo =< y^go^JS^^= A (5 ,^° > 



Rn-l \/Dn-l 



1 ^-Rn-1 <^-Dn-l 

4 -Rn-l -Dn-1 



Finally, 



ED d\ 1 - 1 

res--= <^— 5^.-iA5V'._i>--< 

3+/- 



1 5 Rn-l 5Dn-l 

A „ > 



i? 



n — 1 -^n— 1 



Eres — 



5rfA 1 _, ^ _- 1 5i?n-l , SDn-1 

A -Kn-1 ^ -Kn-1 -t^n-l 



Taking the sum, we obtain 



Uo =< — 5ijn-l A ^V'n-l > 

-n-n-l 
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Lemma 5. (i) The formula 

UQ) = X''<g*JSVA5g> — 

n = ...,—1,0, 1,..., defines meromorphic in Q differential form on T with poles 

at^i,- ■■ ,72Ar and P^'~ ,Pq'~ . 

(a) The symplectic 2-forms defined by the formula 

^n = ^^res^^iQ) 
p 
can be written as 



i ^^ SX 

k=i 



dX. The rest are just 71, ■ ■ ■ , 72JV and Poo and Pq 



Proof, (i) The poles of g* at the branch points (A^, 0) are kiUed by the zeros of 

, The rest are just 71, ■ ■ ■ , 
(ii) By Cauchy's theorem, 

J2 res ^n (<5) + 5^ res e, (g) = 0. 
p k=i 

Near 7^, 



2N 



^;esg^^ 



Therefore, 



SsniQ) = nr^^T^^Hik) + 0(i) = ■ ^;, / A(7fc) + 0(i). 

(A-A(7fc))^ A-A(7fc) 

Note that g* (7/0) = and, using Lemma 3, 

resUQ) = A'^ < gV^yg > A^(7fc)res-^^ = X^{^k)5p{^k) A ^(7^) 
7fe A 7fe A — A(7fcJ A 

We are done. D 

The bracket {•, •j^io is constructed from the symplectic form ujq 

^ ' Sf 5g 6f 6g 






The original AL flow from section 2 can be written 

V'n = {V'n, H}^^, H = N{l2+l2- 2/o " 27o) 

The phase flow is also Hamiltonian 

ijn = {i^n,P}coo, P = N{Io+7o). 
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8. Embedding of the AL system into the function space. First, we 
introduce interpolating trigonometrical polynomials. 

Let f{x) be a smooth 1-periodic complex function and associate to it the 
sequence of interpolating trigonometrical polynomials 

fN{x)= J2 e^^'^'^fNik), iV-2m + l, m = l,2.... 

|A;|<m 

These have the property that /at (a;) = f{x) for x G Sat = {x E E> : x = ^, n = 
0, . . . , A^ — 1}. Let A^AT be a space of such trigonometric polynomials of degree 
m; it is in one-to-one correspondence with the space M^ of N-periodic complex 
sequences: 

fN{x) < > fn, fN{x) e Mn, fn ^ Mn 

if weput e/Ar(f ) = /n, n = 0, . . . , A^ - 1. 
Introduce Hi (ifj) = H^i^ei/j). The region 

B^ = {^^ G Mn •■ NUiiti)) < oo, for all < t < 1} 

is called the "box" of the space A^at; the trigonometric polynomial V'iv belongs 
B^ if and only if |'i/'Ar(ic)| < N, x E Sn- The map iI^n{x) ^^ i/jn allows us to 
define the fiow e^'^s on B^ by the formula 

Evidently B^ is invariant under the fiow ^'^a . Due to the natural embedding 
M-N '^ -M., the dynamics e s also can be defined by the formula e s -^ = 
6**^3 i/jn for any function i/j E A4 which satisfies the inequality |'i/'(a:)| < N, x E 

We define on A^a^ the volume form 



dvol^ 



1 






where V]\[{'i/j) = D^i^i^) and the functional H^ (i/j) is by definition H^^ei/j). 
The fiow e^-^^ preserves both if 3 (-!/;„, -i/j^) and the volume form 
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Therefore, the volume form dvol^ = e '^^dvol and the functional Ti^ are in- 
variant under the flow e^'^s . The finite measure dn^ (i/JjI/j) on B^ with density 

1 



e 2 



■Hi 



^^Uvol"" 



^-^Kw 



NWnW 



y\ i5ip{x) A 5ijj{x) 



xeSf 



is also invariant under the flow e*'^^ . 

To ensure proper analytic control* we need to introduce the function h{x) 
(fig. 5) 



X 

-5^ 



-K -K+1 



K-1 K 



FIG. 5 

and we define the invariant functional X^(^) = h{NTC^)h{—N^T-C^). On the 
box B^ the probability measure dn^i^i/jjifj) 

where S^ is a normalisating factor**. 

Due to the embedding A^at C A4, the measure dfi^i^ifj^i/j) can be transfered 
to the whole of A4. We define on A4 the probability measure duKii^ji') with 
the density 

1 11- 

■^^XK{^)e~ ^'^'' dvol = ^:^XA-(V')e~2'^5 — - A i5ijj{x) A 5ip{x), 



^K 



^K 



OO! 



xes 



where Xk(V') = h{7ii)h{Hs) and S^ is a normalisating factor. 



*Similar cut-off was introduced in [MCV2]. 

**The precise definition of S^ and 'Bk will be given in Lemma 8. 
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In order to prove invariance of the measure d/iK, we use the method of weak 
solutions introduced by McKean, [MC]. The measure d^^ on the initial data 
defines the measure dM^ in the space of paths C {[0,T] ^ M) for T > 0. All 
information about the measure dfx^ and the fiow 6**^3 is encoded now into the 
measure dM^ . The proof of invariance takes three steps. 

Step 1. The measures dfx^ converge to dfiK weakly in H^, 1 < s < |, as 
N ^ oo. 

Step 2. For any fixed K, the family dM^ , N = 1,2... is tight and converges 
to measure dM^- The stationary measure dM^ for fixed t has marginal distri- 
bution dfiK- 
Step 3. The measure dMx is supported on the solutions of the NSL fiow. 

All three steps will be completed in subsequent sections. From the invariance 
of the measure d/iK, it is easy to infer the invariance of the the desired Gibbs' 
state with the density ^e~2'^^dvol. 

9. Convergence of AL Gibbs' state to the NLS Gibbs' state. In the 

previous section, we introduced the family of probability measures dfi^i^i/jji/j) on 
B^ C M with the densities* 

We will show that d^^ -^ d^K, weakly in H^, l<s<|, asA^^oo, where 
dfixii^^i') has the density 

XKWe'^^'^'^^dvol = XK^e-^^'^"^^^ /\idtlj{x) A #(a;). 

oo! ' \ 

xes 

The cut-off in K can be removed easily and d^K converges to d/i. 

We split the integral H^ into a quadratic part G^ (i/j) and a nonlinear part 
R^ (i/j). Using the explicit expression for the integrals /'s, we have 



- /_^(^n+l^n-l + V'n+l'i/'n-l) " m^ni^n+l + -^n^n+l) + 6|^nP 

It is easy to check that 



N-l 
n=0 



*To simplify notation, we omit the normalization factors ^ for a moment. 
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Introduce 

N-l 

|2 



n=0 

and also 

where -0^ = ^ifj^^). With the new notations, 

Now we are ready to prove 

Lemma 6. The family of Gaussian probability measures d'y^^i/j^i/j) with the 
densities 

Nl 



g-,[^^e-W+ivx-W] 1 ^ ^5V(x)A#(x) 



xes, 



N 



converges weakly in H^, ■§ < I' o.^ N ^^ oo, to the probability measure d'y{ip, ip): 



1^ r I ; "|2^| ; |2 

e 2 

xes 



Proof. For a smooth -(/'(s^) = Yl e^'^*^^'i/'(/c) and A^ -^ oo, 

A; 

= A^^e^ ^ 1^ (n + 1/A^) - 2^ (n/N) + ^j {n - 1/N) 1^ + iVe^ ^ |^ (n/A^) 

1 
^ A^^e^ ^ IV'" (n/AT) P + e^ IV' (n/A^) P ^ / IV'"(^)P + mx)\^dx. 



This explains on a formal level, the convergence d'j^ — > ^7. 
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To obtain the full proof, we write all measures in terms of Fourier coefficients. 
First we note that 

1 

\tlj"{x)\^ + \ij{x)\''dx = ^(IGyr^fc^ + lM{k)\^ = ^(72(A;)|V'(A;)|2. 

Using the identity 

N-l 

N 



N-l 



p=0 \k\<m 

for t/j G H^, we have 

|fc|<m 

= J2 cr%{kMM{k)\'. 



k\<m 



Therefore, the measures d'j (ip^ijj) can be written as 



p I fe I < rrt 



\k\<.m 



/\ iS^{k)ASil^{k). 
Similarly, for d'j {t/j , t/j) 



-iE'^ (fe)IV'(fc)l A " ~ 

e fe /\i5ij{k) A5ij{k). 

k 



For any fixed k, 

a%{k) = iV^le^-'^/^ - 2 + e-2-^/^|2 + i ^ iQ^^k^ + i = a\k), 

as A^ ^ oo, whence the convergence of measures on finite- dimentional subspaces 
generated by the Fourier harmonics. 

To prove tightness of the measures 7^, we introduce a "brick" Ba = {i^ : 
|'i/'(/c)| < Ofc, k G Z}, where Ofc > 0. A brick is compact in H^ if and only if 

J2i^ + k^yal<oo. 



*The sign x means up to unessential constant real factor. 
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Let us estimate 



TT ^ 

. -1 9^ 



<7Tl 



27r 



4 



e-^l^(fc)l'l 



5V^(fc) A ^V'lfc) 



\4'{kW<al 



= n 



<7Tl 



1-e" 



4('')"fc 



The inequalities (1 — 6i) x . . . x (1 — 6„) > 1 — 6i 
g-a; ^ p !x~P, X > 0, imply 



bn. bk > 0, and 



7^(5.) >l-J2e- 



^iV^'°^"fc 



|A;|<7 



>l-2Ppl J2 



\k\<m 



1 



Now we obtain an estimate for a%{k) = 4A^"^(cos ^ - 1)^ + 1, |/c| < m. Obvi- 
ously, 

Therefore, 

This, together with 



a'x"^ < (cos a; — 1)^ < b'x"^ for — n < x < n. 

ak^ + 1 <a%{k) <bk^ + l. 



C 



tti 



m + r 



implies 



7™(-Ba)>l 



E 



C>0, 



|/c| + l)"P 



2Pp\ 



Pick any s < f • Then 2s + 1 < 4. Pick any a such that 2s + 1 < a < 4. Then 
the bricks Ba are compact in H^ . Now pick p such that (4 — a)p > 1. Then 
the sum in the last estimate converges. Chose C so large as to make 7^(i?(j) 
arbitrary close to 1. D 

Lemma 7. (i). For any i/j G H^ the functionals 



XW) 



1 



-NX' 



VnW 



'W+JV57^^(V.)] ^ TV 



PAf(V') 



converge to 



< Ve 
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as N ^' CO. 
(n). 

1 > Vn{iI)) > e~^ 
and 

uniformly for all N, provided NHi < K, —N^H^ < K. 
Proof, (i). The nonlinear term 

R^iijn^i'n) =- y^(^n+l^n-l + '^n+l'i/'n-l) |^n |^ 

-6j2\ijn\^ + 12NIo 
can be expressed as 

(A) +^ [(V'n+l - V'n)(V'n - -^n-l) + (^n+1 " V'n)(^n - -^n-l)] IV'" 



n 



|2 



(B) - ^ [(V'n+1 + -^n-l - 2V'n)^n + (^n+1 + -^n-l - 2V'n)^n] IV'^ 

(C) -^ ^ Vn(V'n-l - i^nf + V'n(V'n-l " ^nf 

(D) -X]V'nV'n(^n-l - V'n) + ^^V'n('i/'n-l ' ^n) 

(E) -3 5^ l^nl' - 6 ^ IVnl' + 12iV/o. 

As in the continuous case, if one counts the difference of t/j^s in two neighboring 
points (V'n+i — V'n) and the function t/^n itself of weight 1, then the terms A, B 
and C are isobaric polynomials of degree 6. The term D is isobar ic polynomial 
of degree 5. 

The term A can be reduced to the form 

We transform the terms B and D to a more convinient form also: 

= ^ [(^n+1 - V'n)'i/'n + (-^n+l " V'n)'i/'n] 
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and 

D{lpri,^n) = -^V'n('*/'n-V'n-l)^+V'n(V'n - ^n-lf 
+2j]]|V'n+l-V'nP|V'nP 

The first term in the formula is C with the opposite sign. Finally, 

R{lpn,'tpn) = 

(1) -3^ |^^|4 - 6^^ |^^|2 + 12iV/o 

(2) +^|^,+l-^,_l|2|^„|2 

(3) +'^-Y,\^n-l - ^n?\^n? 

(4) +^^|V'n+l-V'n|'|V'n|' 

(5) + ^ [(-^n+l - ll^n)^n + (V'n+1 " '0n)V'n] 

Q _ 

(6) +- ^ l^n+l - ^n|^ ['i/'nl^n+l " ^n) + ^n(^n+l " ^n)] 

The advantage of such a representation is that the terms (l)-(5) are nonnegative. 
The term (6) is the only term which is not sign-definite. This term vanishes when 
N -^ oo. Finally, as A^ — *> oo 

1 



A^^7^^(V') 


-J 




2\i;' 


NX^{iP) - 


1 


-^' 


PivW -^ 




1. 





(ii). From the definition Hi — N{Io + Iq) = —logD we have Hi {ip) 
— logr'Ar(V'). The inequality NH^ < K implies, that 1 > V]^ >e~^. 
To obtain the lower bound for R^ {^n, V'n), we prove first that 



(7) -3 5^ |^J4 _ g ^ |^^|2 ^ ^2NIo > 2 ^ IV^ 



6 

n I • 
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Indeed, for < a; < 1 

9 S 

- log(l -x)>x + — + —. 
Then, 

-6l0g(l - IV'nP) > 6|^„|2 + 3|Vn|^ + 2|Vn|'. 

This and 12NIq = — Glogf^l — |^nP) imply the inequahty (7). Therefore, the 
term (1) is nonegative. 

We have just one term (6) which is not sign-definite 

Q _ _ 

Using the identity 

and adding to it the two positive terms (3) and (4), we have 
3 V 

2 X] \i^n+l - ^nP {i^ni^n+1 + tpntpn+1 " 'i\'>Pn\^) 
+ 9 5Zl^'^-l -V'nlVnl^ 



1 . 



Therefore, 



|2 



and 



Now we use the constraint —NTi^ < K and the exphcit formula for if 3: 

-H^ii^u, ^n) = Y. I^-+1 - ^-1' - 2|^n|' - 2 l0g(l - |^„|2). 

2 

The estimate — log(l — x)>a; + ^forO<x<l implies 

-21og(l-|^„|2)>2|^,|2 + |^,|4. 
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Therefore, 

Finally, 

K> -N^n^ii;) > ^A^|^(?i + l/A^)-^(n/A^)P + A^-VWA^)^ 

This inequality also provides the estimate for 

max |'0(n/A^)|*. Indeed, 

^J2\'f^{n/N)\^<K, 
implies that there exists n" such that \ijj{n" /N)\ < K^'^. Then, for any n', 

^(n'/AT) = ^ [^(n + 1/A^) - i){n/N)] + i){n" /N) 

n' <.n<,n" 

and, by Schwartz's inquality, 

\^{n/N)\ < J2 mn + l/N) - ^{n/N)\^ + \i^{n /N)\ 



< yjY^ N\ij{n + 1/A^) - V(n/iV)p + \i^{n" /N)\ 

< C2{K). 



Therefore, 



N^n^{t(j) > -4max|V'(n/A^)|2^A^|V'(?i + l/A^)-^(n/A^)|^ 
>-4cl{K)K = -cs{K). 

Also, 

iVJ^(V) = 5^iV-V(n/iV)l' < cl{K). 

Finally, 

W^iij) = -Nl^iilj) + N^n^iij) > -cl{K) - csiK) = -c^K). 
D 



*This is similar to the continuum case when the iJ ^-norm provides a bound in the sup-norm. 
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Lemma 8. (i) The probability measures 

with 






M 

converge weakly in H^, l<s<| as N ^ oo to the probability measure 
with 

M 

(a). The measure d^x converges in H^, l<s<^asK—^ooto the 
probability measure 

with 

M 
Proof, (i). Take M. = H^ x H^ , ^ < s <\ and / bounded and continuous. We 



will 


prove 


that 




(8) 


1 

M 


Vx^ 


1 



M 

as N ^ oo. This implies that S^ -^ Sk as N ^ oo. Then, from (8), we have 



fd^K=^j f^K:^^f-^^''"drii^,^) 



— ^ ^ / fXKe 2^c?7(^,'i/') X ^ = / fd^K X ^ 
as N ^ OO and weak convergence is proved. 
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To begin with, 



fx 



N 



1 -^^''Wa.,Ni 



K ^ < i. e ^" '^'d-i'\^,^) 



^iv(^) 



/ 



^Ar(^) 



1 ) X^e-i^"rf7^+ / /xfe-^^~rf7^(^,^). 



The first term can be estimated using the inequahty of Lemma 8, item (ii): 



<ll/llooe^^^(^) 



M 



Vn{^) 



1 



xUr{^:^) = o{i). 



as N -^ oo. As to the second term, 



(9) 
(10) 



/X^e-^^"rf7^(^,^) = / /xi.e-^^rf7^(^,V) 



A^ / „-^W^ „-^W\ J,, AT 



+ /Xk e 



dr (V', ^) 



+ //(x^-Xi.)e-^^rf7^(V',V'). 



The first integral on the right converges to 



fXxe ^^^d'jiip.ijj). 



The other two integrals vanish as A^ — > oo. Indeed, (9) can be overestimated by 



(11) 



-e^ /X^|W^-WM7^(^,V'). 
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Using the explicit expressions for W and W, we have 



,7V 



(12) 



(13) 
(14) 



-3^ A^|^(n/A^)|^ - 6^ A^3|V'(n/A^)|^ + 12A^%(^) - 2 /" \ipfdx 
+ Y^ N\tjj{n + 1/N) - ^(n - 1/N)\^\ij{n/N)\^ - 4 /" \ij'\^\ij\^dx 



(15) +-^Ar|^(n-l/Ar)-^(n/A^)|2|^(n/Ar)|2-| / \t(;'\^\t(;\^dx 



2 
1 
2 



(16) +i ^iV|V(n + 1/iV) - ^(n/iV)p|^(n/iV)|2 - 1 /" I^Tl^l'rfx 

(17) + J2 N[{ij{n + 1/N) - ^{n/N))i){n/N) + 

+ {^{n + 1/A^) - ^{n/N))ij{n/N)f - l[iP'^ + ^j^j'fdx 



(18) +- ^ A^|^(?i + 1/A^) - ^(n/A^)|2[^(n/A^)(^(n + 1/A^) - ^(n/A^)) + 
+ ^(n/A^)(^(n + 1/A^) - ^(n/A^))]. 



Here, each term contributes nothing when N -^ oo; for example, substituting 
(12) into (11): 



< 



/ 5^iV-i|V'(n/iV)|2-5^ /■ 

J M J n 



i+l/N 



^dx 



n/N 



dr (V, V') 



n+l/AT 
/N 



dr (^, V') 



<A^ 



1/Af 



l/AT 



\ij{x)\^ - \ij{0)\^dx 



dl''{^,^) 



<N / ||^(a:)r-|V'(0)r|rf7''(V'»c^a;. 

J J M 
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Now, by Schwartz, for < |x| < 1/A^, we have 



2 

N 



liv(x)|2-iv(o)nrf7^ </ {mx)\' - mo)\T d^ 

M J JM 



4 j.,W 



M 



2j,,iV 



-2 |^(0)r-|^(x)|1V(0)|^rf7 

'M 



The stationary Gaussian measure d''y^ has the spectral representation 



|fc|<? 



where ijj are complex isotropic Gaussian variables E^\ijj(k)\'^ = a^ (k). Now 

ki,k2 
k3,k4 



The terms in the last sum contributing into the expectation E^ are 

ki = /c2 ki = /c4 ki = k2 = ks = /C4. 

/ca = A;4 /J3 = h 

ki i^k-i kxi^ ks 

Now 

fci/fcs 

+ Yl e'-'^^'^-'-^E^mki)\'E^mk2)\' + J2E''\'P(ki] 

k\^k2 ki 

So, 

E^|V'(o)r-i?^|V(x)p|V(o)p 



= 5Z [l-e'"*"^'^-"^^]^'^|V'(A:i)|'£;'^|^(A;2 



ki^k 



2 



= 5] |l-e2--('=i-'=^)|a^2(^^)^-2(^^)_ 



\ki\<ra 
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Using the estimate a/c"^ + 1 < o"Ar(^) of Lemma 7, uniformly in < a; < 1/N , the 
last sum can be overestimated as 






1 



1 



ak\ + 1 a/c| + 1 



0(1), 



as A^ ^ oo. This imphes that the term (12) contributes nothing when A^ -^ oo. 
To estimate the contribution of (13), expand the log in Xq: 

UN^Ioii;) =-6A^MogJJ(l-iV-2|V'(n/A^)P) 

T\T-4: T\T-6 

=6N'>[N-^\ij{n/N)\^ + -—\^/;{n/N)\^ + —-\^/;(n/N)f 



/^ 



iv) 



4! 



-N-^\il){nlN)\^ 



where f{x) = — log(l — x) and the fourth derivative Z*-*") = 3!(1 — x)~^ in the 
remainder term is bounded since \'i/j{n/N)\ < C2{K) and N~'^\i/j{n/N)\'^ < 9 < 1. 
Now for (13), we have 

1 
25:iV-.|,,VA.)|»-2/|if..4iV-3^/<™.|„ViV)|«. 



The part with the the fourth derivative can be overestimated by 

as N ^ oo. The contribution to (11) of first two terms can be overestimated by 



< 2 



<2N 



pn+l/N 

y^N-^\il){n/N)\^ - / \il)\^dx 

M Jn/N 



d'-y 



N 



l/N 



JM 



(x)r-ii/'(o)r rf7^rfx. 



Again, by Schwarz, 



61 j^.N 



iv(x)r-iv(o)i^u7 



M 



|6 T N 



<2 / 1^^(0)1^^ -|V(0)nV(x)|W =0(1), 

IM 
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as A^ ^ oo, uniformly in < x < 1/A^. This implies that the term (13) makes 
no contribution when A^ -^ oo. The estimation of (14-18) is similar. 
(10) can be overestimated by 



(19) 



XK-Xxld'y 



N 



Now, we have 

xk-xk = [Hnh^) - HHi)] h{-N'n^) + [H-NH^) - h{ns)] HHi). 

Therefore, (19) can be overestimated by 

(20) .-. <||/iu I \h{Nn^) - h{n^)\h{-N^n^)d^'' 

(21) +II/IU j \H-N'n^) - hins)\hini)d^' 



.N 



To estimate (20) we note that \'ilj{n/N)\ < C2{K) on the 
support of h{—N^H^) and that 



NH^ = -NlogY[{l - N-^\ij{n/N)\^) 



Therefore, 

(20)<||/|U||/i'| 



\i>{n/N)\<C2{K) 



< 



Wh'W 



+ 



\i>{n/N)\<C2{K) 



Jo 



Wh'W 

oo "' oo 



\^{n/N)\<C2{K) 



(«) 



EV^"'i^(^/^)i 



rf7^. 



The first term vanishes as in the estimate for the term (20) and the second term 
vanishes as in the estimate for the term (21). 

To estimate the term (21), note that for any e > 0, 



.A^ 



r' Ul^lloo >C) > l-e, 
for all A^ provided C is sufficiently large. Then, 

(21) < 2e+ / \h{-N^n^) - h{ns)\h{ni)d^ 



N 
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in which the integral vanishes for A^ ^ oo as for (20). 
(iv). For any bounded continuous / 



fXKWe--^'^^^^d-f{ilj,i(;) ^ / /e-i^Wrf7(V,V^) 



as i^ -^ oo, by the bounded convergence theorem. This imphes S^ ^ S as 
K -^ oo. Now the last statement is proved. D 

10. Tightness of measures (iiVf^. We assume that the initial data -i/; at (x) has 

distribution dfi^. The solutions of initial value problem ijj(x,t) = 6**^3 ^n{x), 
(x, t) G § X [0, T], are realisations of the complex random field with distribution 
dM^ . To prove tightness of measures we need 

Lemma 9. (Kolmogorov — Centzov). [Ku]. Let the family of complex continuous 
random fields 4>{x,t), {x,t) G S x [0,T], with distributions dm^ satisfy 

(1) i?^i</>(x,t)r<c 

and 

(2) E''\(t){xi,ti)-(t){x2M)V <C{\x^-X2\''' + \ti-t2r'), 

where a^^ +^2'^ < ^ ^'^^ 7 > 0. Then the family of measures dm^ is tight with 
respect to the weak topology of C [(/) : E> x [0, T] -^ C^) 

Since realisations of the field are continuously differentiable in the spatial 
variable x, we consider the random fields i^'{x, t). The main result of this section 
is 

Lemma 10. For any n = 1, 2, . . . ,.■ 

(3) E^\i;'{x,t)\'^<c{n,K) 

(4) E^\ij'{xu ti) - ij'{x2, t2)P" < c(n, K) (\xi - xsP + |ti - tsp/') 

where E^ is expectation with respected to the measure dM^ . 

This, together with Lemma 9, implies tightness of the measures dM'^. Tight- 
ness of the family dM^ can be proved along the same lines. The distribution 
of the limiting stationary random field we denote by dMx- It follows from the 
previous section, that the measure dMx has marginal distributions dfiK- In 
order to establish the Holder continuity of the spatial derivatives of the random 
field we need another 
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Lemma 11. (Kolmogorov-Centzov) . [Ku]. Let (/)(x,t), (x, t) G § x [0,T], he a 
complex random field. Assume that there exist positive constants 7, C, cti and 
a2 with a^ + a^ < 1 satisfying 



E\(l){xi,ti) - (/)(X2, t2)r < C {\xi - X2r + \h - h 



C12'' 



Then the random field i{j{x, t) has a continuous modification. Moreover cq = 
(1 — a^^ — a^^)/7, if /3i and (32 are positive numbers less then aico or a2Co 
respectively, then there exists a positive random constant C with EC^ < 00 such 
that 

\(j)ixi,ti) - (j){x2.t2)\ < C {\xi - X2f' + \tl - t2f') . 

By Fatou's lemma, passing to the limit A^ ^ 00 in the inequality (4), we have 

EkW{xiM) - ^'(X2,t2)|'" < C{n) (|xi -X2r + \ti - ^2!"/') 

Now Lemma 11 implies that the measure dM'j^ is supported on paths satisfying 

\tlj'{xi, ti) - V'(a;2, t2)\ < (\xi - X2\^/^- + \ti - t2|'/^-) . 



One can prove, using methods of [V], that the Holder exponents 1/2— and 1/4- 
are optimal. 

Proof of Lemma 10. To obtain (3), observe that by the invariance of d^^ 



N\j.ti 



E^;!,W{x,t)r = E'^W{xM 



\2n 






by Lemma 7 



Due to Lemma 8, S^ -^ E^ > 0, as A^ ^ 00. Using the Gaussian character of 
the measure ^7"^ the last integral can be overestimated by 



<ci(K,n) 



|2 T .TV 



C2{K,n) 



E k"^~N\k) 



\k\<r 



wm'dY 

LM 

Now the estimate ak'^ + 1 < (7%{k) (see proof of Lemma 6) implies (3). 
To prove (4), write the differential equation for the flow on A^at as in 



difjlx) 



1 



1 



N^{ijj{x + —)+iJj{x - — ) - 2^|J{x)) 



N 



N' 



1 



1 



+\ij{x)Wip{x - — ) + i^{x + — )) 



N- 



N' 
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where if^^x) = "(pN^x^t) G M.n^ x G Sat; equivalently, 



e^-^a [t(;{m)]{x) = e''"^^ [t(j{*)] (x) 
ft 



(5) 
where 



—lie 
'o 



i(t— s)Ajv 



«)r(V'(« + ^,«) + ^(«-^. 



(x)(is, 



AT 



e'i^^[tl;{.)]{x)= ^ e^^'^^'e^'^^^^^^ik), /^{k,N) = N^{u^ + uj'^ -2) 



\k\<7 



and 



\k\<m 



Since (5) holds for all x G Sat with left and right being polynomials, then (5) 
holds for all x eS. 

First, we derive the estimate for spatial increments. Much as before. 



E^\tlj'{x + h,t)-ij'{x,t)\^'' 
= E^\ij'{x + h)-ij'{x)\^^ 

<C2(K,n)[£;^^|^'(/i)-V'(0)lT- 
Using the spectral representation of the Gaussian measure 7^ in the form 



|A;|<m 



we obtain 






The estimate je*^'^^ — 1| < ex, < \x\ < ^ impli 



les 



aA;4 + 1 ^ ^^ ^ aA;4 + 1 

G<\k\<h~^/2 h-^/2<\k\ 



<C3 E 



< 



CAh" Yl 1+C4 Yl k^-""^^- 

0<|fc|<ft,-i/2 h-'^/2<\k\ 
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Finally, 

(6) E^\i;\x + h,t)- ij\x, t)|2- < C6(i^, n)/i". 

To obtain the estimate for temporal increments we use (5) and write* 

E^\i;'{x,t + h)- ^'(x, t)|2- = E^lij'ix, h) - ^'(x, 0)|2- 

(7) <CTE^\e'^^''^'{x)-^'{x)f" 

h 



(8) +cr E^ 



Here 



J{h—s)Aff 



1 



1 



'2'V(«+^,s)+^(.-^,s) 



{x)ds 



N 



2n 



(7) < eg E^N\e'''^^ij'{x)-ij'{x)\^'' = C9 [E^N\e'^^^ij'{x)-ij'{x)\^]" . 
Now, from the spectral respresentation, 

|A;|<7Ti 

(9) = J2 {i-coshN''{u'' + u-''-2))\n^ea-''{k) 



|fc|<m 



(10) 



+ J2 sin^ hN^iiv'' + UJ-'' -2)4Tr^k^a]^'^ {k). 



|fc|<m 



First, we note that for |x| < tt 

(11) ciqx'^ < 1 — cosx < ciix'^ 
Therefore, for \k\ < m, we have 

(12) ci2k'^ <-N'^{cj'' + u-'' -2)<ci3k^. 

To estimate (8), we assume that \k\ < h~^''^a^ where a is a suitably chosen 
constant. Then, by (12), 

\hN^{u^ +u-^ -2)\ <7r. 



*„itA 



e'^^N^I^x) = e^^^N [^(,)] (a;) 
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Using (11) and (12), we find 

|fc|<m, 
\k\<h-^/'^a 



|fc|<m 



< 



ci&h^ E k^ + cn E -^ < cis/i 



1/2 



|fc|</i-i/2Q, |A;|>/i-i/2a 

The estimate of (10) is even simpler: 

(10) <c,, Y. ^+-- E ^ 

|A;|<m |fc|<m 

|A;|</i-i/2 |fc|>/i-i/2 



<C2 



l/l^ E /c" + C22 Y^ p < C23/1 



1/2 



Finally, for (7) we have the estimate 
(13) E^\e''^^iP'{x)-iP'{x)\^'^ <C2^{K,n)W''^. 

To estimate (8), observe that 



E 



J{h-s)AT 



1 



1 



\^|;{.,s)\\^l;{. + -,s)+^|;{.--,s)) 



N 



N- 



{x)ds 



N 



2n 



< /i^-i / E^ 



J{h — s)AN 



1 



1 



m;s)\\i;{. + -,s) + i;{.~-,s)) 



N 



N' 



2n 



ds 



<h 



2n-l I Tj^N 
K 



E 



AsA.r 



IV'(-,«)P(V'(- + ^,«) + V'(--^,«)) 



N 



N 
2n 



ds. 



We will show that, for any s > 0, 



(14) E 

Then, 

(15) 



itAjv 



1 



1 



|2(^(. + -,.)+^(. --,.)) 



N 



2n 



< C25{K,n). 



^N 



2n 



(8)=Ej^| ...ditto... r" <C25/l^". 
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The estimates (13) and (15) produce 

E^liP'ix, t + h)- ij'ix, t)|2- < C26{K, n)h''/^ 

This and (6) produce (4). 

To prove (14) first we note, is 'iIj{x) = 'Y^'\p[k)e^'^'^^^ , then 

^n{x) = J2 i^N{k)e'^'''\ with ^N{k) = Yl ^(^)- 

|fc|<Tn n=k, (mod AT) 

To simphfy notation we write t(j = i/jn, then 



1 



1 



x)\^ {ip{x+ — ) +ilj{x- — ) 



A^ 



A^' 



2nk? 



J2 e2v^^(fcl+fc3-fc2)x2 ^^g _^^^k^^^^k^^^^k 



\ki\<m 



and 



^itAjv 



m*r{^i*+j;^)+^i 



1 



N 
2nk 

'n 



(x) 



N 



J2 e^tA(N,p)^2.^pik,+ks-k,).2 cqs —±2mp{ki + k^ - k2)^{ki)^{k2)'4^{k:i] 
\ki\<m 



where p{k) = k (mod N) and \p{k)\ < m. 
Now, 



drr.e 



itAr 



1 



1 



iV'(«)rU'(- + i^)+^(--i^ 



A^' 



A^- 



(x) 



N 



^ V^ git[A(7V,p)-A(Ar,p')]g27rix(p(A;i+fc3-A:2)-p(fc4+fc6-fc5))2cos^^2cOS^^ 

^ N N 

\ki\<m 

X An'^p{ki + k3- k2)p'{ki + ks - /C5)V'(^l)^(^2)'i/'(fe)V'(/S4)V'(/j5)V'(/^6) 

Estimating 



E^N ... ditto . . 



l2n 



\E^n\ . . .ditto. 



|2n| 



y^ TT ^it[A{N,p^)-A{N,p'^)]^27.ix{p{k{-ki-k2)-v{k{+kl-ki))2 ^^g ^^^3 q ^^g 27r/c^ 



J = l 



A^ 
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X E^NY[47r^p{k{ + ki- ki)p'{ki + kl- kl)^{ki)^{ki)ij{ki) 

X ^{ki)tl;{ki)^{ki)\ 

n 

^E^7- n "^^'iPiK + ki-kmPiK + ki-ki) 
i=i \ki\<m 

Using the inequalities 

\p{ki + ks-k2)\< \k^ + k2- ks\ < (1 + |A;i|)(l + |A;2|)(1 + |A;3|), 

and 

(1 + |A;|)2<2(1 + |A;|2), 

we overestimate the last sum by 

\pi\<m 

The worst term 

J2E,NplmPi)\' pLlV'(3n)|2. 

can be overestimated by 

.•.<E^iv|V''(a:)px3n<C(?i,A^) 

due to a/c^ + 1 < a^(A;). D 

11. Identification of measure dM.K- This section shows that the measure 
dM-K is supported on the solutions of the NLS flow. 
Multiplying both parts of the original equation 

i^* = _^" + 211/^1 V 
by the test function /(x,t), (x, t) G S x (0, T] and integrating produces 

/ dx dt [it(jf - t(jf" + |V'|^2V'/] = 0. 

J-oo 



AN ADDITIONAL GIBBS' STATE 
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Lemma 12. In the statistical ensemble dAix 



Emi 



1 p+(x> 
J -oo 



dxdt[ii;f -i;f" + \ij\^2^f] 



0. 



Proof. The equation for the AL flow on A^tv is 



iij'{x,t) = -A7v^(x,t) + \ij{x,t)\'^{ij{x+ — ,t) + ^(x- — ,t)) 



N' 



N' 



where An'i(j{x) = N"^ (V'(x + ;^) + ^("^ " 'k) ~ '^'^(^)) ^^^ ^ ^ ^a^- ^^ multiply 
by the test function f{x,t) and integrate over the variables t and x. 



1 r+oo 



J-oo 



dx dtdN{x)[ii{j{x, t)f* — i^{x, t)A]\[f 



1 



1 



+ |^(a;, t)\'[ijix+—,t) + ij{x-—,t) 



N' 



N' 



in which 5Nix) = 4? Yl ^(^ ~ ^o)- Now 



xo6Sj> 



Emi 



J Jdxdt[tijr-^f" + \i^\^2i^f] 



N- 



lini i?A^iv|... ditto ...I 



^x 



lim -Eyvi^l / / dx dt iijjf*{l — 6n{x)) — / / dx dt ip{d^^ — 5NAN)f 



+ 




dx dt 



IV'|'2V' 



1 



1 



< lim cE 



N- 



M\ 



+ lim cE 




2 

dx dti/j f* {1 — Sn 



/P 



JV- 



xi 




dxdtip{6'^^ — (JatAat)/ 



+ lim ci? 



AT- 



M{ 




dxdt 



|V|'2V' - IVI' ( i^ix + h + i^{x - i-) ) dN{x) 



N' 



N' 



f 
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The first term overestimated by 



.■.<\r\ooEj^N\ I / dxdt\tp{x,t)-tp{k/N,t)\dx 

<\r\ooT j E^nU \i;{x,t)-ij{k/N,t)\dx 
<|rUT2 I dxEMN\^{x,t)-^{k/N,t)\^ 

.l/JV 

=\r\ooT''N \ £;^^|V(x)-V(o)pdx 

■J 

0<x<l/Ar 

as A^ ^ oo due to the stochastic continuity of the random field "0. Estimates for 
the remaining terms can be obtained along the same lines. D 

Due to the results of [MCVl, B2] the initial data ilj{x,t) determines the fiow 
?/'(•, t) for all t. In probabilistic language this means that '(/'(•, t) is measurable 
with respect to the field generated by ?/'(•, 0). 



[AL] 
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